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We have investigated the phase transition properties of 
classical linear sigma model. The fields were kept in contact 
with a heat bath for sufficiently long time such that fields 
are equilibrated at the temperature of the heat bath. It was 
shown that the sigma model fields undergoes phase transition, 
but the transition temperature depend crucially on the lattice 
spacing. In the continuum limit, the transition temperature 
tends to zero or at least to a very low value. 
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The possibility of forming disoriented chiral conden- 
sate (DCC) in relativistic heavy ion collisions has gen- 
erated considerable research activities in recent years. 
Anslem and Ryskin |Q first considered the possibility of 
producing classical pion fields in heavy ion collisions. The 
idea became popular after Rajagopal and Wilczek g] pro- 
posed the quench scenario and Bjorken et al ||] proposed 
the Baked-Alaska model for producing DCC. Rajagopal 
and Wilczek argued that for a second order chiral phase 
transition, the chiral condensate can become temporarily 
disoriented in the nonequilibrium conditions encountered 
in heavy ion collisions. As the temperature drops below 
T c , the chiral symmetry begins to break by developing 
domains in which the chiral field is misaligned from its 
true vacuum value. The misaligned condensate has the 
same quark content and quantum numbers as do pions 
and essentially constitute a classical pion field. The sys- 
tem will finally relaxes to the true vacuum and in the 
process can emit coherent pions. Since the disoriented 
domains have well defined isospin orientation, the associ- 
ated pions can exhibit novel centauro-like fluctuations 
of neutral and charged pions. 

Most dynamical studies of DCC have been based on 
the linear sigma model PH^l in which the chiral degrees 
of freedom are described by the real 0(4) field <& = (a, tt), 
with the Lagrangian, 



(i) 



where A is a positive coupling constant and /„■ is the pion 
decay constant. In the vacuum the symmetry is sponta- 
neously broken and the sigma field acquires a non zero 
vacuum expectation value < 0|cr|0 >= / w . Symmetry 
is restored by a second order phase transition. Transi- 
tion temperature can be calculated in mean field level 
as, T c 2 = 12 fl/{N + 2), with N=4 0. In this model pi- 
ons are massless. To be realistic, one introduces a small 



symmetry breaking term H a to the lagrangian and pion 
become massive. The parameters of the model can be 
fixed by defining the meson masses and pion decay con- 
stant. 

It is well known that classical fields suffers from di- 
vergence problem. Thermodynamic of classical fields are 
only defined if an ultraviolet cut off is imposed on the 
momentum. In lattice simulations, finite lattice spac- 
ing provides the cut off. Some thermodynamic quanti- 
ties such as the diffusion rate of the topological charge 
Jl8| , p~9[ | in non-abelian gauge theories, or the Lyapunov 
exponent p(| which is equivalent to the damping rate 
of soft thermal excitations are found to be insensitive to 
the cut off. However, it is not known, whether the critical 
temperature of phase transition the linear sigma model 
depend on the momentum cut off or not. With the possi- 
bility of experiments being performed at RHIC in search 
of Quark-gluon Plasma, where one expects to find signal 
of DCC, this investigation becomes important. 

Aim of the present letter is to show that critical tem- 
perature of linear sigma model do indeed depend on 
the lattice spacing or the momentum cut off and in the 
continuum limit, the phase transition temperature ap- 
proaches zero. To see phase transition we have used 
the old idea of finite temperature field theory, that is if 
you attach the fields with a heat bath, their equilibrium 
configuration will corresponds to that of the heat bath. 
We thus solve the equation of motion of sigma model 
fields in contact with heat bath which we represent as 
a white noise source. To be consistent with fluctuation- 
dissipation theorem, we also include a friction term in the 
equation of motion. We thus propose to study following 
Langevin equation for 0(4) fields, 



[a + v d t + \(<f> 2 - tf)]<S> = Hn a + ( 



(2) 



where rj is the friction. The heat bath £ is represented by 
a white noise source with zero average and correlation as 
demanded by fluctuation-dissipation theorem, 



< C(t,x,y,z) >= (3a) 
< Ca(h,xi,yi, z 1 )Cb(t2,X2,U2, z 2 ) > d 4 x = 2Ti]S ab (3b) 



Recently it was shown that in </> 4 theory hard modes 
can be integrated out on the two loop level leading to 
dissipation and noise in the quasi classical limit for the 
propagation of the long wavelength fields This also 
justify our approach. 



If the fields are allowed to evolve in contact with the 
heat bath at temperature T for a long time, they will be 
equilibrated at that temperature. Above a certain critical 
temperature, the fields will undergo symmetry restoring 
phase transition, if the model allows for such a transition. 
For sigma model, the condensate value of the sigma field 
provide a convenient way to probe the phase transition. 
In the symmetric phase, the condensate is zero, while in 
the symmetry broken phase it is non-zero. 

In the present paper we will consider linear sigma 
model without the symmetry breaking term. It is well 
known that there will not be an exact phase transition 
if the symmetry breaking term is included. We solve the 
eq|2] on a 32 3 lattice with periodic boundary condition. 
Solving eqn.|^ require initial conditions (<fi and </>). We 
distribute the initial fields according to a random Gaus- 
sian with 



< cr > - < a >=< tt? > 



(4a) 
(4b) 

< m > l = v z /4f(r) (4c) 
< cr >= < 7r 4 >= (4d) 
(4e) 



< cr >= f(r)f„ 

< TTi > = 
2_ „.2 , 



< <T 2 > = < 7T > 2 = V 2 



The interpolation function 

f(r) = [1 + exp(r - r^fT)}- 1 



(5) 



separates the central region from the rest of the system. 
We have used r$ = 11a where a is the lattice spacing and 
r=0.5 fm. The initial field configuration corresponds to 
zero temperature, < cr >~ f v , < tt >^ 0. Initial zero 
temperature fields will thermalise to the temperature of 
the heat bath if kept in contact for sufficient time. The 
other parameter of the model is the friction (rj). In the 
present paper, we use rj = r/ n + r/ a and for r/ n and rj a use 
values as calculated by Rischke jl3j but its precise value 
is not of importance here, as we are looking for fields at 
equilibrium. Friction determines the rate of approach to 
equilibrium. This aspect of equilibration was verified. 

We define volume averaged sigma condensate < a > 
as the order parameter, 



< a >= 1/V J 



d 3 x<7 



(6) 



In fig. 1, we have shown the equilibrium value of the 
order parameter as a function of the temperature. Re- 
sults for 4 different lattice spacing, a=0.5, 1.0, 2.0 and 2.5 
are shown. For all the lattice spacings, the order parame- 
ter decreases from ~ / w at very low temperature to exact 
zero at some critical temperature (T c ) and then remain so 
beyond the critical temperature. Phase transition in the 
model is evident. However it is also evident that phase 
transition temperature (the temperature at which a con- 
densate vanishes) depend strongly on the lattice spacing 
used. It goes to smaller and smaller values as the lattice 



spacing is reduced. It seems that in the continuum limit, 
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1. Order parameter as a function of temperature for 



four different lattice spacing, a=0.5, 1.0, 2.0 and 2.5. 
outer ones are for higher lattice spacing 



The 



The present result indicate that classical sigma model 
fields donot show correct equilibrium behaviour. Phase 
transition temperature depend strongly on the lattice 
spacing, and in the continuum limit, it goes to zero or 
atleast to a very low temperature. Strong dependence 
of the the transition temperature on lattice spacings in- 
dicate that the simulation studies of disoriented chiral 
condensate using classical sigma model fields are highly 
questionable. The results obtained from those simula- 
tions can not be believed to represent physical systems. 
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